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I. INTRODUCTION 

The detailed study of top quark pair production close to threshold constitutes a major 
part of the top quark physics program at the International Linear Collider (ILC) during 
its first phase running at lower energies. While the top quark mass (in a threshold mass 
scheme |l|, |2|, |3() can be measured with a precision of around 100 MeV from the cm. energy 
where the cross section line-shape rises, accurate measurements of the top quark total width, 
the strong coupling and (in the case of a light Higgs) the top Yukawa coupling can be 
gained from the normalization and the form of the line-shape j^|. Electroweak effects play 
a crucial role for theoretical predictions of top pair threshold production. Apart from the 
ti production process, mediated either by virtual photon or Z boson exchange, one also 
needs to account for initial state QED beam effects and the finite top quark lifetime in the 
leading order approximation. On the other hand, due to the fact that, close to threshold, the 
top quarks are produced with small relative velocities, u<l, singularities proportional to 
(a s /v) n and (a s mt>) n arise from virtual and real gluon radiation in 0{a™) QCD diagrams 
that need to be summed up to all orders in a s . These singularities arise from ratios of 
the kinematical scales m t ("hard"), pt ~ m t v ~ m t a s ("soft") and E t ~ m t v 2 ~ m t a 2 
( "ultrasoft" ) that are governing the nonrelativistic ti dynamics, where pt and E t are the 
top three-momentum and kinetic energy, respectively. Since the Standard Model top quark 
width T t ~ T t (t — > bW) ~ 1.5 GeV is numerically comparable to the average top kinetic 
energy, it serves as an infrared cutoff and protects the nonrelativistic ti dynamics from 
nonperturbative QCD effects. 

In recent years most work in the literature on top pair threshold production was focused 
on determining higher order QCD corrections using nonrelativistic QCD (NRQCD). 1 In the 
fixed-order approach, where the singularities oc (a s /v) n are systematically summed, next- 
to-next-to-leading order (NNLO) corrections to the total cross section are fully known |l| 
and even some NNNLO analyses have become available recently (see e.g. Refs. p,|6j]). In 
the renormalization group improved approach both singularities oc (a s /v) n and oc (a s lnv) n 
are systematically summed. Here, the full set of NNLL order corrections is known except for 
the full NNLL running of the dominant S-wave ti production current p, |g, |9( . At present 
the theoretical QCD normalization uncertainties of the NNLO fixed-order prediction for the 
total ti cross section close to threshold is estimated to be of the order 20% [l|, while for the 
(yet incomplete) NNLL order renormalization group improved prediction it is 6% (lOj. 

As far as electroweak corrections are concerned no complete determination of all effects 
beyond the leading order level has yet been achieved for the total cross section. At leading 
order the three electroweak effects mentioned above can be treated independently: The top 
pair production is described by (e + e~)(ti) effective NRQCD operators with Wilson coeffi- 
cients that account for the propagation of the virtual photon and the Z boson mediating the 
ti production process. The QED beam effects consisting of the machine-dependent beam 
energy spread, beamstrahlung and initial state radiation can be conveniently implemented 
through a luminosity spectrum that is convolved into the theoretical predictions without 
any initial state QED corrections ^, |ll|, |l2j. The top quark width, on the other hand, is 
incorporated into NRQCD by an imaginary matching condition of a bilinear top quark mass 
operator and can be effectively accounted for by shifting the cm. energy into the complex 



1 We use the term NRQCD to refer to a generic low-energy effective theory which describes nonrelativistic 
tt pairs and bound state effects and not for a theory valid only for scales mt > (J, > mtv. 



plane, y/s — > y/s + iT t [l3|, [lj|. In Ref. [15j a systematic procedure for the theoretical 
treatment of the effects from the finite top quark lifetime for the total cross section be- 
yond leading logarithmic (LL) order in the renormalization group improved approach was 
presented. It was shown that the effects of the top decay can be systematically incorpo- 
rated when NRQCD is matched to the Standard Model by including the absorptive parts 
in electroweak matching conditions that are related to the top quark decay A systematic 
expansion scheme can be achieved by using the power counting |l5l . llfij | 

v ~ a s ~ g ~ g' (1) 

for a s and the SU(2) and U(l) gauge couplings g and g' . The parametric counting scheme 
in Eq. (pQ) is motivated by the close numeric relation between the average top quark kinetic 
energy and the top quark width, E t ~ m t a 2 s ~ T t ~ m t a. It shows that one-loop electroweak 
corrections to the tt production process can contribute at NNLL order. In particular, it was 
shown in Ref. |l5( that, due to gauge cancellations, there are no electroweak matching 
corrections to the NRQCD potential operators and heavy- quark-gluon interactions up to 
NNLL order. However, nontrivial (real and imaginary) electroweak matching conditions 
exist for the dominant (e + e~)(tt) operators as well as for the top bilinear operators due to 
higher order corrections to the top quark width and the relativistic time dilatation effects. 
In Ref. |l5j the imaginary parts of these electroweak matching conditions were analyzed in 
detail. It was found that they lead to corrections of the line-shape form between 2% and 
10% and can shift the IS peak position by up to 50 MeV. 

In this work we analyze the real parts of the NNLL electroweak matching conditions 
of the dominant S-wave (e + e~)(ti) operators. We exclude pure QED corrections 2 from 
our considerations. A complete treatment of pure QED effects together with a consistent 
computation of NNLL initial state (beam) and final state corrections shall be postponed 
to subsequent work. The corrections considered in this work have been discussed in detail 
before in Ref. [l7| and were partly reanalyzed again later in Ref. |la j. Here we present the 
results of a new independent computation and point out a number of discrepancies to the 
results in these references. 

The program of this paper is as follows. In Sec.|H] we explain the notation relevant for 
the matching conditions of the (e + e~)(tt) operators considered in this work. In Sec. lIHI we 
present our results and discuss the differences to Refs. |l7[ and |l8[ . Section |IV| contains a 
brief numerical analysis, and in SecMwe conclude. 

II. NOTATION 

In this section we set up the notation for the electroweak matching conditions to the 
dominant S-wave ti production and annihilation operators relevant for e + e~ collisions. To 
be specific we use the notations that were introduced in Ref. [l5j in the framework of 
vNRQCD J19l|20j]. The dominant effective theory operators to be used for tt pair production 



2 In this context pure QED corrections correspond to loop diagrams that contain a photon but no other 
Standard Model gauge bosons. 



in an S-wave spin triplet state have the form 

Ov, P = [e l3 e]0^ 11 
O a , p = [e 7 ;75e]C£i 

where 



(2) 



(3) 



and the indices j = 1,2,3 are summed. They give the contribution AC = 
J2 P (CVty) @v,p + Ca{v) Oa, p ) + H.c. to the effective theory Lagrangian where the hermi- 
tian conjugation is referring to the operators only and does, due to unitarity, not act on the 
complex Wilson coefficients [15j- The term v is the dimensionless vNRQCD renormalization 
group parameter, where v — 1 corresponds to the hard matching scale (/x = m t ) and v ~ a s 
indicates the range of scales where the effective theory matrix elements are evaluated. The 
relation between the Wilson coefficients at the matching scale v = 1 and for v < 1 has 
the multiplicative form Cv,a(v) / CV,a(1) = exp(/(m 4 , v)). The function / is fully known at 
NLL order |20L |21 |. At NNLL order only the nonmixing contributions to the anomalous 
dimensions J9j and the subleading mixing effects from the spin-dependent potentials J22J are 
presently known. The matching conditions for the Wilson coefficients at NNLL order (and 
without accounting for pure QED corrections) can be written in the form (i = V, A) 
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The terms Cy° A n refer to the tree level matching conditions originating from virtual photon 
and Z boson exchange and read 
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Qf is the fermion electric charge, s w (c w ) the sine (cosine) of the weak mixing angle and 
a is the electromagnetic coupling. The factors multiplying the tree level contributions 
represent the hard QCD matching conditions which are presently known up to order a 2 . 
The term C VA ' a s in Eq. (J3J) describes the imaginary (absorptive) parts of the matching 
conditions. They arise from cuts in the one-loop electroweak diagrams that are related to 
intermediate states that arise in the top decay and were identified in Ref. |l5| . In predictions 
for the total cross section they account for interference contributions of the amplitude for the 
double resonant signal process e + e" — > ti — >• bW + bW~ with the single resonant amplitudes 
describing the processes e + e~ — > t + bW~ — > bW + bW~ and e + e" —> bW + i —> bW + bW~. 
According to the power counting scheme in Eq. (0) these interference effects contribute at 
NNLL order to the NRQCD matrix elements for the total cross section. In Ref. [l5| it was 



also shown that these NNLL matrix elements contain phase space UV-divergences that lead 
to the NLL running of (e + e~)(e + e~) forward scattering operators. Finally, the terms C™ A 
in Eq. (0J) refer to the real parts of the NNLL (one-loop) electroweak matching corrections, 
which are the main focus of this work. 

Using the optical theorem the contribution of the effective theory operators in Eqs. (J2J) 
to the total ti production cross section close to threshold reads 
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where (k + k' = (y/s, 0) and e = k/|k|) 
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= ~(k + k') 2 (5 lk -e l e k ) 
is the spin- and angular-averaged lepton tensor and (q = (y/s — 2m t , 0)) 
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is the time ordered product of the ti production and annihilation operators O 3 x and O x , 
which describes the nonrelativistic dynamical effects. In the case of nonzero electron or 
positron polarization the formula in Eq. (jHJ) is simply modified by 
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where P T refers to electron/positron polarization along the momentum direction. Since Cy 
and Ca have opposite signs (see Tab.|TJ) the cross section is large for positive e + and negative 
e~ polarization. In the unpolarized case the relative NNLL corrections to the total cross 
section coming from the real parts of the one-loop electroweak matching conditions read 
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where all coefficients can be taken at v — 1. 



III. MATCHING PROCEDURE AND RESULTS 



The NNLL order electroweak matching conditions for the Wilson coefficients of the effec- 
tive theory operators in Eqs. (J2J are obtained from the one-loop electroweak corrections to 
the Standard Model amplitude for e + e~ —* ti for on-shell external top quarks in the limit 
y/s — > 2m t (i.e. v — > 0), where the top quarks are at rest. For simplicity the electron mass 
is neglected except for the self energy corrections, the CKM matrix is considered to be the 
unit matrix and the bottom quark mass is neglected in the jti and Zti vertex corrections 
(including the top quark wave function renormalization). If QED corrections are neglected, 



there are no contributions singular in v that are associated with nontrivial NRQCD ma- 
trix elements and all contributions need to be absorbed directly into the effective theory 
matching conditions. Neglecting pure QED corrections the Standard Model amplitude can 
be written in the form 
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where u)± = |(1 ± 75) and k + k' = 2p = (2m t , 0). The real parts of the NNLL electroweak 
matching conditions then read 
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To our knowledge the coefficients h™ R were first computed in Ref. [17j in Feynman gauge 
using the on-shell scheme for the initial state e + e~ pair, the final state tt pair and the 
intermediate Z boson, and employing the fine structure constant for the electromagnetic 
coupling a -1 = 137.036. We have determined the real parts of the coefficients Hl,r in the 
same renormalization scheme by hand and, independently, using the automated packages 
FeynArts [23[ and FormCalc J24J (except for the Z'j box diagrams). For the presentation 
of our results we follow (with minimal modifications) the conventions and notations used in 
Ref. [17[. Since all components of the coefficients h^^ were presented explicitly in Ref. |l7j 
we discuss in detail only those elements where our results differ. For undefined symbols and 
variables used in the following discussion we refer the reader to Ref. LLj 
The results for the coefficients h'£' R can be cast into the form 
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and 

& = -£&■, /}[ = *» ~ slQ > , (21) 

t 3 j being the third component of the weak isospin of fermion /. 

The terms h^ E R describe the self energy corrections where LL^,, i = AA, ZZ, ZA, are 
the renormalized transverse photon and Z self energy corrections and the photon-Z mixing 
correction. Compared to the expressions given in Ref. [17[ our result does not sum the photon 
and Z boson self energy corrections into the denominator of the respective propagators 
since the resulting higher order corrections are beyond NNLL order and because at the ti 
threshold the intermediate photon and Z boson are far off-shell. The terms h e L R describe 
the corrections to the e + e~ vertex, where F % L R , i = A, Z, are the vertex corrections to the 
left- and right-handed e + e~7 and e + e~Z vertices. The terms h l [ R refer to the corrections to 
the ti vertices, where the coefficients a A,z correspond to the corrections to the fit and Zti 
vertices, respectively. The arguments of the various terms in Eq. (J20JI indicate the virtual 
bosons that are being exchanged in the various triangle diagrams, e.g. a A,z (W, W) refers to 
the vertex diagrams with exactly two internal W lines and a A,z ((pz) to those where the only 
internal boson is the neutral pseudo-Goldstone boson. Note that a A,z (W, 4>w) > a> A ' z {.<l>z,H) 
and a A (Z,H) are identically zero in Feynman gauge. Finally, the terms h h ^ R describe the 
contributions from the WW, ZZ and Z'-f box diagrams. We note that the Z7 box diagrams 
are infrared-finite for the ti pair being at rest. 

Our results agree with those given in Ref. [17[ except for the following points: 3 

1. For the corrections to the left-handed e + e~7 and e + e~Z vertices the coefficients F L ' 
read 

Ft = £ (Qe (PI? P(Q, M z ) - JU(g, M w )^ , (22) 

Ft = f ((Plf p(q,M z ) + (f3l-2tf^)^p( q ,M w ) -^A^Mw)) .(23) 

47T y \ S w / ^S w S w S w J 

We agree with the result for the function p in |l7j], but for the function A we find 
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For A cw , which describes the corrections to the total cross section from the real parts 
of the NNLL order electroweak matching conditions, this leads to an absolute shift by 
+0.012 with respect to the results in Ref. |l7| |. 



We acknowledge that the discrepancies 2) and 3) were confirmed to us by [2 



2. For the corrections to the tt vertices coming from the exchange of the neutral physical 
Higgs boson we find the result 

M 2 
+ B 1 {\P, M H , \M) - 3j^B (P, ±M, \M) \ , 

where the upper component refers to the photon and the lower to the Z vertex, respec- 
tively. Our results are consistent with the ones in Refs. |l8|,l26|; but differ from those 
in Ref. |l7| in the argument of the last B function. Furthermore the counterterm 
contributions read 

<4 = —{Qt{SZ R + SZ L ) + 24 6k) , (27) 
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4 = - (WZ R + fifiZ L ) + 2^4 8k) . (28) 

a \ s w J 

For the left- and right-handed wave function renormalization terms we have the ex- 
pression 

5Z LtR = A LtR + m 2 t (A' R + A' L - 2C") , (29) 

where C = Jp-(p 2 = m 2 ) . For the function C we find 
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which is consistent with Ref. [18|], but differs from Ref. [17[, where the B\ function was 
obtained instead of the first B§ function. For A ew these changes lead to an absolute 
shift by +0.076 with respect to the results in Ref. (l2( for m# = 130 GeV. 

3. For the corrections from the Z7 box diagrams we find the result 

hp = ;£(-/$)(& - ffiQeQW* 1 , (31) 

hV = ^IWr ~ PDQeQaF** > (32) 

where (u = -^f-) 

F Zl = -^—(B (lP,^M,0) + B ( 1 1 P,lM,M z )-2B (P,0,M z )) 

(Jj X 

+ (B {P, 0, M z ) - B Q {1, 0, \M)) + (u- l)M 2 C (l, \P, 0, §M, M z ) . (33) 



The expression for F Zl agrees with the one in Ref. J22J, but differs from the one in 
Ref. |l7[ by an additional minus sign in front of the first term. We note that the same 
error is also contained in the analysis of Ref. |l8j . For A ew the changes lead to an 
absolute shift by +0.004 with respect to the results in Ref. |l7| . 

8 



At this point we also take the opportunity to point out a number of discrepancies we 
find to Ref . |l8[ . (For the definition of the various undefined variables we refer the reader 
to Ref. Ilal.) For the WW box contribution we find 
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This differs from |18| (in the published version) by the factor r multiplying the function 
/i and the global minus sign in the definition of j\. For the contribution of the Z'j box 
diagrams we find the expression 

F^(r) = + 21nr + lnl-(-l-^(ln(l-l)-^- I l-/ 4+ g-l) / 3.(36) 

For the tt vertex correction a ZH , which corresponds to the triangle diagram containing a Z 
boson and a physical neutral Higgs boson, we find the expression 
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which differs from the one in [l^J by a minus sign in front of the function /g 
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IV. NUMERICAL ANALYSIS 



In this section we give a brief numerical discussion of the real parts of the NNLL elec- 
troweak matching conditions obtained in this work. In Tab. |1] the numerical values for Cy°™ 
and C™ A are displayed for various values for the top and the Higgs masses, a" 1 = 137.036 
and 



M, 



w 



m P 



m,, 



m r 



= 80.425 GeV, 


M z 


= 91.1876 GeV, 


c 2 = 


= M^/Ml , 


0.511 MeV, 


«v 


= 0.106 GeV, 


m T 


= 1.78 GeV 


0.005 GeV, 


m d -- 


= 0.005 GeV, 


m s 


= 0.10 GeV 


1.3 GeV, 


m b -- 


= 4.2 GeV 







(38) 



for the gauge boson, the lepton and the quark masses. Note that the finite electron mass is 
applied in the calculation only for the self energy corrections. 

The vector coefficients dominate for the tree level as well as for the one-loop coefficients. 
The one-loop corrections show a significant Higgs mass dependence and vary in the vector 
(axial- vector) case between 8.6% (—1.2%) and 5.6% (—2.6%) for Higgs masses between 
115 GeV and 1 TeV and m t = 172.5 GeV. 



mt (GeV) 


170 


175 


m H (GeV) 


115 
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200 


1000 
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C born( 10 -7 GeV -2) 


-5.429 








-5.123 








C born( 10 -7 GeV - 2 ) 


1.260 








1.184 








C^(10- 8 GeV- 2 ) 


-4.562 


-4.073 


-3.702 


-3.060 


-4.460 


-3.951 


-3.566 


-2.890 


C^ w (10- 9 GeV- 2 ) 


-1.416 


-2.286 


-2.797 


-3.025 


-1.408 


-2.335 


-2.904 


-3.260 


pew /rtborn 


0.0840 


0.0750 


0.0682 


0.0564 


0.0871 


0.0771 


0.0696 


0.0564 


/oew /^~iborn 


-0.0112 


-0.0181 


-0.0222 


-0.0240 


-0.0119 


-0.0197 


-0.0245 


-0.0275 



TABLE I: Numerical values for the tree level and real one-loop electroweak matching conditions 
and Cy* A , respectively, and the C v w A /Cy° A Tl ratios for various values for the top and the Higgs 



/-"vborn 

masses, a~ l = 137.036 and the values given in Eq. (|38|) . The coefficients Cy ™ do not depend on 



the Higgs mass. 
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FIG. 1: The correction A cw as a function of run for mt = 172.5 GeV and the values given in 
Eq. (|38|) for a -1 = 137.036 (dashed curve) and in the scheme for the electromagnetic coupling 
defined in Eq. (|3*9*|l (solid curve). 



A more transparent view on the impact of the real electroweak one-loop corrections on the 
predictions of the total ti threshold cross section can be gained by considering the quantity 
A ew denned in Eq. (|12|). In Fig. ^ the dashed line represents A ew as a function of the Higgs 
mass for m t = 172.5 GeV and adopting the previous choices for the other parameters. For 
m H = (115, 150, 200, 1000) GeV we have A cw = (0.161, 0.142, 0.128, 0.104). The Higgs mass 
dependence is rather strong for small tuh and drops quickly close to the decoupling limit for 
increasing m#. The strong Higgs mass dependence for small Higgs masses originates from the 
fact that for a light Higgs boson the dominant effects of the virtual Higgs exchange between 
the ti pair could be described in NRQCD by a Yukawa potential that leads to a singularity 



oc m t /mH for nix — > (see e.g. Refs. [18|,l26j,|28|])- In an approach where the Higgs exchange 
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would be accounted for in NRQCD through a Yukawa potential the electroweak one-loop 
matching conditions would need to be modified to account for the contributions caused by 
the Yukawa potential in the NRQCD matrix elements. If the Higgs boson is indeed very 
close to the lower LEP bound this would be a viable alternative to our approach where all 
virtual electroweak effects are encoded in the NRQCD Wilson coefficients. 

The solid line in Fig. Q finally, shows A cw for the same choice of parameters but adopting 
an MS definition for the QED coupling at the scale \i = m t that accounts for the leading- 
logarithmic vacuum polarization effects due to the three charged leptons and the quarks 
below the top quark scale, 
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In this scheme a nf 8 (/z = mt) = 1/125.926 for m t = 172.5 GeV and the real NNLL elec- 
troweak matching conditions are modified, 

C"r = C5S, - Cfc ^^ ( £ Ql In (§) + E NM In K) ) . (40) 

\i=e,/i,r ^ l ' i=u,d,c,s,b ^ l / 

This scheme is preferred to the one where the fine structure constant is used, since it 
leads to a substantially reduced dependence on the light fermion masses in A ew . This 
is because most of the fermionic vacuum polarization effects are cancelled by the mod- 
ification in Eq. (pfU|) . Moreover in this scheme A ew is by 17% smaller than in the one 
with the fine structure constant. For m# = (115, 150, 200, 1000) GeV we now have 
A ew = (-0.001,-0.021,-0.037,-0.063). The results show that concerning the real parts 
of the NNLL order electroweak matching corrections (up to the QED effects neglected in 
this work) the effects from virtual Higgs exchange are comparable to the other electroweak 
effects. 



V. CONCLUSION 

We have determined the real parts of the NNLL order (one-loop) electroweak matching 
conditions for the Wilson coefficients of the dominant NRQCD operators describing ti pro- 
duction close to threshold in e + e~ annihilation. The results include the contributions of 
all electroweak one-loop effects that are integrated out when NRQCD is matched to the 
Standard Model except for pure QED corrections. The results are an important ingredient 
for a complete treatment of electroweak effects at NNLL order for top quark pair production 
at a future Linear Collider. We have pointed out a number of discrepancies with respect to 
earlier work. 
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